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We calculate radiative corrections to the Lamb shift of order a s (Za) 5 m and 
radiative corrections to hyperfine splitting of order a 3 (Za)Ep generated by 
the diagrams with insertions of radiative photons and electron polarization 
loops in the graphs with two external photons. We also obtain the radiative- 
recoil correction to hyperfine splitting in muonium generated by the diagrams 
with the t polarization loop. 



Nonrecoil corrections of order a 3 (Za) 5 m to the Lamb shift and corrections of order 
a 3 (Za)Ep to hyperfine splitting are generated by three-loop radiative insertions in the 
skeleton diagram in Fig. 1. Respective corrections of lower orders in a generated by one- 
and two- loop radiative insertions are already well known (see, e.g., review [1]). The crucial 
observation, which greatly facilitates further calculations, is that the scattering approxima- 
tion is adequate for calculation of all corrections of order a n (Za) 5 m and a n {Za)Ep (see, 
e.g., a detailed proof in [2]). One may easily understand the physical reasons which lead to 
this conclusion. Consider the matrix elements of the skeleton diagram in Fig. 1 with the on 
shell external electron lines calculated between the free electron spinors, and multiplied by 
the square of the Schrodinger-Couiomb wave function at the origin. They are described by 
the infrared divergent integral 
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in the case of the Lamb shift, and by the infrared divergent integral 1 



8Za dk 
hip 



Tin 3 



in the case of hyperfine splitting. In these integrals k is the dimensionless momentum of the 
exchanged photons measured in the units of the electron mass. 



FIG. 1. Skeleton two-photon diagram 

Let us consider radiative insertions in the skeleton two-photon diagram in Fig. 1. Account of 
these corrections effectively leads to insertion of an additional factor L(k) in the divergent 
integrals above, and while this factor has at most a logarithmic asymptotic behavior at 
large momenta and does not spoil the ultraviolet convergence of the integrals, in the low- 
momentum region it behaves as L{k) ~ k 2 (again up to logarithmic factors), and improves 
the low-frequency behavior of the integrand. However, the integral for the Lamb shift 
is sometimes still divergent after inclusion of the radiative corrections because the two- 
photon-exchange diagram, even with radiative corrections, contains a contribution of the 
previous order in Za. This spurious contribution should be removed by subtracting the 
leading low-momentum term from L(k)/k A . The result of such subtraction is a convergent 
integral, where the low integration momenta (of atomic order mZa) in the exchange loops 
are suppressed, and the effective loop integration momenta are of order m. Then it is clear 
that small virtuality of the external electron lines would lead to an additional suppression 
of the matrix element under consideration, and it is sufficient to consider the diagrams only 
with on-mass-shell external momenta for calculation of the contributions to the energy shifts. 
As an additional bonus of this approach one does not need to worry about the ultraviolet 
divergence of the one-loop radiative corrections. The subtraction automatically eliminates 
any ultraviolet divergent terms and the result is both ultraviolet and infrared finite. 

Below we consider contributions to the Lamb shift and hyperfine splitting generated 
by radiative insertions in the skeleton diagram in Fig. 1. We also obtain radiative-recoil 
correction to hyperfine splitting generated by the r polarization loop. 



We define the Fermi energy Ep as 

EF = fz*a^(l + a»)(^y c hR 00 , (2) 

where m is the electron mass, M is the muon mass, m r is the reduced mass, a is the fine structure 
constant, c is the velocity of light, h is the Planck constant, is the Rydberg constant, a M is the 
muon anomalous magnetic moment, and Z is the nucleus charge in terms of the electron charge 
(Z = 1 for hydrogen and muonium). 
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II. CORRECTIONS OF ORDER a 3 {ZafM TO LAMB SHIFT AND OF ORDER 

a 3 {Za)E F TO HYPERFINE SPLITTING 



A. Diagrams with Three One-loop Electron Vacuum Polarizations 



1. Lamb Shift 



Each polarization loop in the diagrams in Fig. 2 corresponds to insertion of the vacuum 
polarization operator (a/n)k 2 Ii e in the Lamb shift skeleton integral in eq.(l), where 



he 



f 

Jo 



dv 



v 2 (l-v 2 /3) 
4 + k 2 (l - v 2 ) 



(4) 



Inserting also the multiplicity factor 4 we obtain an analytic expression for the contribution 
to the Lamb shift generated by the diagrams in Fig. 2 in the form 



SEP = — m/ dkk 2 l\ e . 

7r 4 n 3 V m ) Jo 



Calculating the integral numerically we obtain 



SE^ = - 0. 021 458 (1) 



a 3 (Za) 5 (m 



7r 2 n 3 V m 



— m, 



(5) 



(6) 



or 



SE^ = -0.002 16 kHz 



(7) 



for the IS* level in hydrogen. 



+ 



FIG. 2. Three one-loop polarizations 



2. Hyperfine Splitting 

We obtain the expression for the radiative correction to hyperfine splitting generated by 
the diagrams in Fig. 2 inserting the polarization loops in the skeleton integral in eq.(3) 
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After numerical calculations we have 



6E { jp FS = 2. 568 3 (4) 



a 3 (Za) 



(9) 



or 



5E { h ] fs = °- 003 29 kHz 



(10) 



for the ground state in muonium. 



B. Diagrams with Two-Loop and One-loop Electron Vacuum Polarizations 

1. Lamb Shift 

The integral for the diagrams in Fig. 3 is obtained from the skeleton integral in eq.(l) 
by insertion of the one-loop vacuum polarization (a/7r)k 2 I le , and the two-loop vacuum 
polarization (a/n) 2 k 2 I 2e (see, e.g., [3,4]) 
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where the dilogarithm Li 2 (x) is defined as Li 2 (z) = — J* 2 dt\n(l — t)/t. 

Inserting in the skeleton integral in eq.(l) also the multiplicity factor 6 we obtain an 
analytic expression for the contribution to the Lamb shift generated by the diagrams with 
the one- and two-loop polarization blocks in Fig. 3 



5E? = 



96a 3 {Za) 5 /m r \ 3 



7r 4 n 3 



m 



roc 

m / dkli e l 2e - 
Jo 



(12) 



After numerical calculations we obtain 



5E 



(2) 



0.390 152 (7) *l§& ,„. 



(13) 



or 



5E { ? = -0.039 21 kHz 



(14) 



for the IS level in hydrogen. 
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FIG. 3. One- and two- loop polarizations 



2. Hyperfine Splitting 

In the case of hyperfine splitting we obtain the expression for the energy shift generated 
by the diagrams in Fig. 3 with the help of the skeleton integral in eq. (3) 



48a 3 (Za) 
After numerical calculations we have 



je^Ej-^t (15) 

7r 4 n d Jo 



8E$ FS = 3. 559 9 (2) E F , (16) 

7T 



or 



SE^ls = 0.004 56 kHz (17) 
for the ground state in muonium. 



C. Diagrams with Three-Loop Electron Vacuum Polarization 

1. Lamb Shift 

For calculation of the correction generated by the diagrams in Fig. 4 we need the three- 
loop vacuum polarization operator {a/n) 3 k 2 I^, e . This operator in QED and QCD was con- 
sidered in a series of papers [5-9]. As a result seven leading terms both in the low- and 
high-momentum asymptotic expansions in the powers of the momentum were calculated 
analytically. Some of the coefficients were presented in [7,9] only in the MS scheme and 
only for the case of QCD. We adjusted these results for the case of the momentum renor- 
malization scheme used in QED, and constructed an interpolation which approximates the 
three-loop polarization operator for all Euclidean momenta. 

The skeleton integral in eq.(l) remains infrared divergent even after insertion of the 
three-loop vacuum polarization since i3 e (0) ^ 0. This linear infrared divergence is effectively 
cut off at the characteristic atomic scale mZa if we restore finite virtualities of the external 
electron lines. As was already mentioned in the Introduction, such infrared divergence lowers 
the power of the factor Za, and respective would be divergent contribution turns out to be 
of order a 3 (Za) 4 . This correction was calculated in [10], and we will not discuss it here. We 
carry out the subtraction of the leading low-frequency asymptote of the polarization operator 
insertion, which corresponds to the subtraction of the leading low-frequency asymtote in the 
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integrand for the contribution to the energy shift h e (k) = he(k) — he(®), and insert the 
subtracted expression in the formula for the Lamb shift in eq.(l). We also insert an additional 
factor 2 in order to take into account possible insertions of the polarization operator in both 
photon lines. Then the contribution to the energy shift has the form 



Calculating the integral numerically we obtain 



SE^ = 1.015 88 (5) - 77 [ — ) m, (19) 



a 3 (Za) 5 / m r \ 3 



n 2 n 3 v m 



or 

,(3) 



5E£> = 0.102 10 kHz (20) 
for the IS level in hydrogen. 




FIG. 4. Three-loop polarizations 



2. Hyperfine Splitting 

In the case of hyperfine splitting there is no problem of infrared divergence for the 
radiative correction generated by the three-loop polarization insertions in Fig. 4. This 
correction is given by the integral 

7r 4 rr jo 

which arises after insertion of the doubled three-loop polarization operator in the skeleton 
integral in eq.(3). 

After numerical calculations we obtain 

5E$ FS = 1.647 9 (5) E F , (22) 



or 

n(3) 



5E ( h' fs = 0.002 11 kHz (23) 
for the ground state in muonium. 
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D. Diagrams with One-Loop Electron Factor and Two One-loop Electron Vacuum 

Polarizations 



1. Lamb Shift 

An analytic expression for the correction of order a 3 (Za) 5 generated by the gauge in- 
variant set of diagrams in Fig. 5 can be obtained from the skeleton integral in eq.(l) in 
the same way as the other corrections above. But this approach requires knowledge of a 
new element, namely, the gauge invariant electron factor L L {k) in Fig. 6 which describes all 
possible insertions of the radiative photon in the electron line with two external photons. 
An explicit expression for the electron factor was obtained in different forms in [11-14] (we 
use the expression from [14]) 



r 1 !, ,2 1 ^ 



Ink — In 



2k 



(24) 
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$(k) = k 



dz 



o 1 — k 2 z 2 



In 



1 + k 2 z(l -z) 
k 2 z 



(25) 



Inserting in the skeleton integral in eq.(l) the electron factor (a/ir)k 2 LL(k), one-loop 
polarization operator squared and the multiplicity factor 3 we obtain the radiative correction 
in the form 



A8a 3 (Za) 5 (m^ 3 



7r 4 n 3 



m 



dkk 2 L L (k)I 2 le . 



(26) 



It is easy to check explicitly that this integral is both ultraviolet and infrared finite. The 
infrared finiteness nicely correlates with the physical understanding that for the diagrams 
in Fig. 5 there is no correction of lower order a 2 (Za) 4 generated at the atomic scale. 
After numerical calculations we obtain 



SE^ = 0.0773 (4) ( ^ 



a 3 (Za) 5 (m r 



7r 2 n 3 V m 



m, 



(27) 



or 



5E ( l ] = 0.007 77 (4) kHz 



(28) 



for the IS* level in hydrogen. 
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FIG. 5. One-loop electron factor and two one-loop polarizations 




FIG. 6. One-loop electron factor 
2. Hyperfine Splitting 

We calculate the contribution to hyperfine splitting generated by the diagrams in Fig. 5 
using an explicit expression for the electron factor like in the case of the Lamb shift above. 
This is a different electron factor which corresponds to a different spin projection. It was 
obtained in [15] and has the form 

3 4 , (2 1 In A; 2 1 
LnFs{k) = -- 2 -- 2 \nk 1T ^ + - 



. , 2 VF+4, VW+l + k 

In A; ; In 



k 



(29) 



9v& 2 + 4 Vk 2 + A + k 4 , Vk 2 + 4 + k l, n . 4_, J , 

H 7~, In ; ; hi , h —r$(k) - T^$(fc). 

2 k 3 ^fWTl-k WWTl y/W+l-k ^k y ' k 3 

Inserting the electron factor (a/ir)k 2 LHFs{k) together with the one-loop polarization 
operator squared and the multiplicity factor 3 in the skeleton integral in eq.(3) we obtain 
the radiative correction in the form 

SE&s = ^y^ EF r dkk 4 L HFS (k)ll (30) 
After numerical calculations we obtain 

8E$ FS = - 3.487 2 (2) E F , (31) 

or 

8E { h FS = -0.004 47 kHz (32) 
for the ground state in muonium. 
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E. Diagrams with One-Loop Electron Factor and Two-Loop Electron Vacuum 

Polarization 



1. Lamb Shift 

An integral representation for the correction generated by the diagrams in Fig. 7 is 
obtained from the skeleton integral in eq.(l) in the standard way 



32a 3 (Za) 5 fm 



7r 4 n 3 



— ml dkL L (k)I 2e . 
m J Jo 



(33) 



Calculating this integral numerically we obtain 



5Ef = 2.191 3 (4) 



a 3 (Za) 5 fm r \ 3 



it 2 n 3 \ m 



m, 



(34) 



or 



SE^ = 0.220 24 (4) kHz 



(35) 



for the IS level in hydrogen. 



FIG. 7. One- loop electron factor and two- loop polarization 



2. Hyperfine Splitting 

We obtain the radiative correction to hyperfine splitting generated by the diagrams in 
Fig. 7 inserting the electron factor (a/7r)k 2 LfjFs{k) together with the two-loop polarization 
operator and the multiplicity factor 2 in the skeleton integral in eq.(3) 



lQa 3 (Za) 
7i 4 n 3 

After numerical calculations we obtain 



SEjPps = ^y^ EF r dkk 2 L HFS (k)he. (36) 
7r*n° Jo 



6eP fs = -4.680 9 (1) - E F , (37) 
or 

SeP fs = -0.006 00 kHz (38) 
for the ground state in muonium. 



a 3 (Za) 
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F. Diagrams with One-Loop Polarization Insertions in the Electron Factor and in 

the External Photon 



1. Lamb Shift 

The contribution to the Lamb shift generated by the diagrams in Fig. 8 is similar to 
the contribution in eq.(33), the only difference is that now we consider a radiatively cor- 
rected electron factor in Fig. 9 and a one-loop polarization insertion in the external photon. 
Insertions in the skeleton integral in eq.(l) lead to the expression 



32a 3 {Za) 5 fm 



7r 4 n 3 



^) m f°° dkL^\k)I le , 
m J Jo 



(39) 



where the parametric representation for the electron factor with one-loop polarization inser- 
tion in Fig. 9 has the form [13] 



tf'(^/ i v'';-y' L t (M). 

Jo 1 — 1> 



(40) 



where A 2 = 4/(1 — v 2 ), and L L (k, A) is the one-loop electron factor in Fig. 6 for a massive 
photon with mass A. An explicit representation for this electron factor was obtained in [13] 
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d(x, A) = x 2 + A 2 (l - x), 



(42) 
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a 2 (x, y, A) = d(x, A) + k 2 y(l - y), 
Calculating the integral in eq.(39) numerically we obtain 



SEf } = 0.037 36 (1) 



a 3 (Za) 5 (m 



7r 2 n 3 v m 



or 



5Ef } = 0.003 75 kHz 



(43) 
(44) 



for the IS level in hydrogen. 



FIG. 8. One-loop polarization insertions in the electron factor and external photon 



7 T 



4^ +2 4^ + ^ 



FIG. 9. One-loop polarization insertions in the electron factor 



2. Hyperfine Splitting 

We obtain the radiative correction to hyperfine splitting generated by the diagrams in 
Fig. 8 inserting the radiatively corrected electron factor (a j 'ir)k 2 L^'p S (k) in Fig. 9 together 
with the one-loop polarization operator and the multiplicity factor 2 in the skeleton integral 
in eq.(3) 

SE f Fg = 1 ^§^E F / o ~ dkk^(k)I le . (45) 

The parametric representation for the electron factor L H ' FS (k) has the form [16] 

L$p s (k) = f 1 dv v2{1 ~ v2 J 3) L HFS (k, A) , (46) 
Jo 1 — v z 

where A 2 = 4/(1 — v 2 ), and L}jFs{k, A) is the one-loop electron factor in Fig. 6 for a massive 
photon with mass A. An explicit representation for this electron factor was obtained in [16] 

r u = i/', f x H ( A{X-x,y) k 2 B(\;x,y) ^ 

HFS[ ' } 2 Jo X Jo V \k 2 y(l - y) + x 2 + A 2 (l - x) [k 2 y(l - y) + x 2 + A 2 (l - x)} 2 J 

(47) 
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where 



and 



A 2 (l - x) 

A(X;x,y) = a (x,y) + <n(x, y) ^ - ^ _ , (48) 
B(X; x , y) = b„(x, y) + bl (x, y) - + „) ( _^A_ J , (49 ) 



a (x, y) = (l-xf-x- 2— + -(1 - -)y\ (50) 

XXX 



ai (x, y) = (-- 3(1 -x))y + {- 2 --- 2)y 2 , (51) 

\x J X z X 



b (x, y) = x(l - l)y + (-- + 1 + x)y 2 + (A _ 1 _ 3 )y 3 + -y 4 , (52) 
2 a; ar x x 



r 2 m 8 2 

h(x, y) = (- - 1 - 2x + -)y 2 + (-- + - + 4 - 2x)y 3 + (1 - -)y 4 , (53) 

37 ^ *2T 3? 37 



4 - x 
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M^!/) = -r(i-i)f (54) 

x z 



After numerical calculations we obtain 



5eP fs = -0.533 3 (5) - ^ E F , (55) 

or 

= -0.000 68 kHz (56) 

for the ground state in muonium. 

G. Diagrams with Two One-Loop Polarization Insertions in the Electron Factor 

1. Lamb Shift 

The contribution to the Lamb shift generated by the diagrams in Fig. 10 is similar to the 
correction generated by the one-loop polarization insertion in the electron factor calculated 
in [13]. The explicit expression for this correction 

w = _itow (! m* m r<ft 4» \ k) - 4»'(o) 

7r 4 rr \ m J Jo k 
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differs from the respective expression in [13] only due to the difference between the elec- 

(2 1) 

tron factor with one one-loop polarization insertion L L ' (k) in eq.(40) (see Fig. 9) and the 
electron factor with two one-loop polarization insertions L L ' (k) in Fig. 11. 
The photon line with two one-loop polarization insertions has the form 
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where X\ = 4/(1 - v 2 ) and X 2 = 4/(1 - v 2 ). 

Then the electron factor with two one-loop polarization insertions in Fig. 11 can be 
written as 
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Calculating the integral for the energy shift in eq.(57) numerically we obtain 
5EP = -0.012 610 (3) a3{Za)5 ^ 



TT 2 n 3 



m 



or 



SEP 



-0.001 27 kHz 



for the IS level in hydrogen. 



FIG. 10. One-loop polarization insertions in the electron factor 
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FIG. 11. Two one-loop polarization insertions in the electron factor 
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2. Hyperfine Splitting 



The contribution to hyperfine splitting generated by the diagrams in Fig. 10 is similar to 
the correction generated by the one-loop polarization insertion in the electron factor which 
was calculated in [16]. The explicit expression for this correction has the form 



8a 3 (Za) 
7i A n 3 



E, 
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where (compare eq.(59)) 
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After numerical calculations we obtain 
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7T^ 



or 



6E$ FS = -0.000 40 kHz 
for the ground state in muonium. 
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H. Diagrams with Two-Loop Polarization Insertion in the Electron Factor 

1. Lamb Shift 



The contribution to the Lamb shift generated by the diagrams in Fig. 12 is similar to 
the correction generated by the one-loop polarization insertion in the electron factor which 
was calculated in [13]. The explicit expression for this correction 



se(s) = _lQa\Zaf 
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differs from the respective expression in [13] only due to the difference between the electron 

(2 1) 

factor with one-loop polarization insertion L L ' (k) in eq.(40) (see Fig. 9) and the electron 



(3 2) 

factor with the two-loop polarization insertion L L ' (k) in Fig. 13 
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where A 2 = 4/(1 — v 2 ), and the electron factor with a massive photon L L (k,\) is written 
explicitly in eq.(41). 

A convenient expression for the subtracted massive electron factor L L (k,\) — L L (0,\) 
was obtained in [13], and using those old formulae we immediately obtain 



^ = -0.245 71 (7)5^^^, m 



n 2 n 3 \ m 

or 



SE ( P = -0.024 70 kHz (69) 

for the IS level in hydrogen. 



FIG. 12. Two-loop polarization insertions in the electron factor 



i i = i r + 2 J — H — r + -r — 
+ 2 r ? r + 4 J — H — r + 2 -r — 



FIG. 13. Two-loop polarization insertions in the electron factor 



2. Hyperfine Splitting 



The contribution to hyperfine splitting generated by the diagrams in Fig. 12 is similar to 
the correction generated by the one-loop polarization insertion in the electron factor which 
was calculated in [16]. The explicit expression for this correction has the form 
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+ 



3 1 - v 2 
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and LnFs{k, A) is the electron factor with a massive photon from eq.(47). 
After numerical calculations we obtain 



5E§ ] FS = -0.123 9 (6) 



o?{Za) 



E h 



7T Z 



(72) 



or 



6E$ FS = -0.000 16 kHz 



(73) 



for the ground state in muonium. 



III. r POLARIZATION CONTRIBUTION 



The one-loop r-lepton polarization contribution to hyperfine splitting generated the dia- 
grams in Fig. 14 may be calculated exactly. Again the scattering approximation is sufficient 
for calculation of this correction (see, e.g., [17]). First time the r-lepton contribution was 
estimated in [18]. At that moment this correction was of purely academic interest, and a 
crude step-function model for the one-loop polarization spectral function was used in [18]. 
Due to a spectacular experimental progress during the last two decades, now we need a more 
accurate result for the r-lepton contribution to hyperfine splitting. 

The general expression for this correction has the form (compare, e.g., [2]) 
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(74) 



(75) 



is the one-loop r-lepton vacuum polarization, the dimensionless parameter \i is given by the 
expression \i = m e /(2m M ), and the Fermi energy E F unlike the expression in eq.(2) does not 
include the factor 1 + a M . The expression in Eq. (74) may be obtained from the integral for 



the skeleton graphs with two exchanged photons by the substitution 1/k 2 
additional factor 2 has the combinatorial origin. 



2I lT , where the 



+ 2 



2: 



FIG. 14. r lepton polarization contribution to hyperfine splitting 



After the Wick rotation, transition to the four-dimensional spherical coordinates, and 
to the dimensionless integration momenta measured in the units of the electron mass the 
expression in eq.(74) acquires the form 
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8E T = 



a(Za) m e 
n 2 m„ 



roo 

/ dkk 


\-( 


h 


.fik \ 



-K^-*) + 5(i^-T-5 



1 , v 2 (l -v 2 /3) 
o *4 + k 2 (l- v 2 ) 

mi \ / 



a(Za) m e 

06t Cjp. 



7T" 



m. 



This is a finite integral which can be calculated numerically with arbitrary accuracy. Nu- 
merically we obtain 



5e T = 0.019 190 6.... 



(77) 



One can also obtain an analytic expression for leading terms in the expansion of the r- 
lepton polarization contribution to the hyperfine splitting over the small parameters m M /m T , 
m e /m^, and m e /m T . Let us describe briefly calculation of the leading terms in this expan- 
sion. First, we write the dimensionless contribution to the energy splitting as a sum of two 
terms 



5 €l = 8 J™ dkk — ^l + fi 2 k 2 -fikj-^ (pky/l + /i 2 k 2 - fi 2 k 2 



(78) 



Se 2 = 8 



f 



dkk 



^ + k 2 (l-v 2 ) 



L 



i , v 2 (l-v 2 /3) 
dv- 



4m% 



f + k 2 (l-V 2 ) 



(79) 



which correspond to the two first and two last terms in the square brackets in the integrand 
in eq.(76), respectively. The integral Se2 is proportional to m 2 /m 2 , and is too small to be 
of any interest for us here. The integral 5ei, as we will see, is proportional to a much larger 
parameter m^/m 2 , and gives a leading contribution to 5e T . To calculate it we once again 
rescale the integration momentum q = fik 



i , v 2 (l-v 2 /3) 
dv - — -2 . 

fe) +? 2 (1-^ 2 ) 
(80) 



To extract the leading terms in the asymptotic expansion of this integral we introduce 
an auxiliary parameter a which satisfies the inequality 1 <C a <C rrir/m^. We use the 
parameter a to separate the momentum integration into two regions, a region of small 
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momenta < q < a, and a region of large momenta a < q < oo. In the region of small 
momenta we use the low momentum expansion of the polarization operator and obtain 



S e < = 8 r dq 
Jo 



1 + q 2 - q) - 



Q 



_4_ / V 2 

15 \m T , 



(81) 



/6, 6 , 1\ / m u 
- In a + - n2 + - — ^ 

V5 5 2/ U T 



In the region of large momenta g>lwe use the large momentum expansion of the skeleton 
integrand and obtain 



5e>= f°° dq (P\ C dv-^ 
J a \2qJ Jo f Hi 



v 2 (l -v 2 /3) 



fe) +fa-v 2 ) 



(82) 



6 m T 6 , 6, n 77 
- In ma In2H 

5 m„ 5 5 50 



m. 



For the intermediate momenta q ~ cr both approximations for the integrand are valid simul- 
taneously, so in the sum of the low-momenta and high-momenta integrals all a-dependent 
terms cancel and we obtain a cr-independent result 



5e\ = 5e^ + de^ 



5 m„ 25/ \m T 



0.019 185. 



(83) 



The leading terms in the asymptotic expression for the r-lepton contribution were also 
estimated in [18]. We disagree with the both terms obtained in [18]. However, numerically 
for the real parameters of the r-lepton, the difference between the result in [18] and in eq.(83) 
is only about 4 x 10~ 3 . 

Comparing the approximate result in eq.(83) with the result of the numerical calculation 
of the integral in eq.(76) we see that their difference is about 5 x 10~ 6 . Due to overall 
smallness of the correction under consideration, this means that the analytic expression in 
eq.(83) is sufficient for all phenomenological purposes. 

Finally, the r polarization contribution to the hyperfine splitting may be approximated 
by the expression 



8E T = 5e 1 



a(Za) 



7T^ 



m,. 



6 , m 7 
-In — 

5 m. 



+ 



51 \ a(Za) m e m 



25 



7T 



mi 



Hip, 



or numerically 



5E T 



0.002 2 kHz 



(84) 



(85) 



for the ground state in muonium. 
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IV. DISCUSSION OF RESULTS 



In this paper we calculated a series of nonrecoil corrections of order a 3 (Za) 5 m to the 
Lamb shift, and a series of nonrecoil corrections of order a 3 (Za)Ep to hyperfme splitting 
generated by the diagrams in Figs. 2, 3, 4, 5, 7, 8, 10, 12. Collecting all contributions to 
the Lamb shift in eq.(6), eq.(13), eq.(18), eq.(27), eq.(34), eq.(43), eq.(60), and eq.(68) we 
obtain 

SE» = 2.651 9 (6)23M (IM 3 m , (86) 

or 

5E t ° t = 0.266 53 (6) kHz (87) 

for the IS level in hydrogen. 

Collecting all contributions to hyperfme splitting in eq.(9), eq.(16), eq.(21), eq.(31), 
eq.(36), eq.(55), eq.(64), and eq.(72) we obtain 

5Ef FS = - 1.358 (1) E F , (88) 

or 

5Ef FS = -0.001 74 kHz (89) 

for the ground state in muonium. 

Both the corrections to the Lamb shift and hyperfme could be easily estimated before 
the actual calculation is carried out. They are suppressed by an additional factor a /n in 
comparison with the corrections of the lower order in a. In the case of the Lamb shift this 
means that corrections of order a 3 (Za) 5 m should be as large as 1 kHz for the IS level in 
hydrogen. Corrections of this magnitude are phenomenologically relevant at the current 
level of the experimental and theoretical accuracy (see, e.g. [1]). We expect that the largest 
contribution will be generated by the gauge invariant set of diagrams with insertions of three 
radiative photons in the electron line in the skeleton diagrams in Fig. 1. Work on calculation 
of the contribution of these diagrams as well as of all other remaining corrections of order 
a 3 (Za) 5 m to the Lamb shift, and corrections of order a 3 (Za)Ep to hyperfme splitting is 
now in progress, and we hope to report on its results in the near future. 

We also obtained above the r lepton polarization contribution to the hyperfme splitting 

SE^fUn^ + ^)^1'^E F , (90) 
\5 m M 25/ n 2 m 2 



which numerically gives 



SE T = 0.0022 kHz (91) 



for the ground state in muonium. 

The magnitude of this contribution is comparable to a number of other new corrections, 
obtained recently, for example to some nonlogarithmic three-loop radiative-recoil corrections 
[19], and to the contributions due the two-loop hadron polarizations in [20]. 
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